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HANKEL DETERMINANT APPROACH TO GENERALIZED VOROB’EV-YABLONSKI 

POLYNOMIALS AND THEIR ROOTS 

FERENC BALOGH, MARCO BERTOLA, AND THOMAS BOTHNER 


Abstract. Generalized Vorob’ev-Yablonski polynomials have been introduced by Clarkson and Mansfield 
in their study of rational solutions of the second Painleve hierarchy. We present new Hankel determinant 
identities for the squares of these special polynomials in terms of Schur polynomials. As an application of 
the identities, we analyze the roots of generalized Vorob’ev-Yablonski polynomials and provide formulae for 
the boundary curves of the highly regular patterns observed numerically in [8]. 


1. Introduction and statement of results 
Let u = u{x; a) denote a solution of the second Painleve equation 

Uxx = XU 2u^ H- a, x € C. 


(1.1) 


It is known that for special values of the parameter a C C the equation admits rational solutions. In fact 
Vorob’ev and Yablonski [23, 24] showed that ioi a = n G Z, the equation (1.1) has a unique rational solution 
of the form 


u{x] n) = — In 
da: 


Qn-ijx) 

Qn{x) 


n G Z>i; 


u{x\ 0) = 0, 


u{x;—n) =—u{x; n), n G Z>i, (1-2) 


which is constructed in terms of the Vorob’ev-Yablonski polynomials {Qn{x)}n>o- These special polynomials 
can be defined via a differential-difference equation 


Qn+i{x)Qn-i{x) = xQlix) - 4:(^Q'^{x)Qn{x) - {Q',^{x)y^ , nG Z> 1 , xGC, 
where Qq{x) = 1, Qi(a:) = x, or equivalently [18] in determinantal form: with qk{x) = 0 for A: < 0, 

4 


e„(*) = nSdet 


k=l 


2^k\ 


in-l 

qn- 2 i+j{.x) , nGZ>i- ^qk{x)w = exp 

J l,j—0 

A;=0 


— -VJ -\-wx 

o 


(1.3) 


(1.4) 


For our purposes, it will prove useful to rewrite (1.4) in terms of Schur polynomials. In general (cf. [20]), 
the Schur polynomial s\ G C[t] in the variable t = .. .),tj G C associated to the partition A = 

(Ai, A 2 ,..., A^(a)) with Z 3 Xj > Aj+i > 0 is determined by the Jacobi-Trudi determinant, 

SA(t) = det [/iA,-j+fc(t)]^|^l^ . (1.5) 

Here, hfc(t) for k G Z>o is defined by the generating series 


^hk{t)z’" = exp 


and /ifc(t) =0, fc < 0. 


( 1 . 6 ) 


k=0 


va=i 


Remark 1.1 (Homogeneity). From (1.6) it follows immediately that hk{t) is a weighted-homogeneous func¬ 
tion, 

/ifc(t) = e'^hk (e“4i, ■ ■ ■) , eG C\{0}, 
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and hence also 

i(\) 

sx{t) = e^^^sx{c~^ti,e~‘^t2,e~^h,..) , |A| = ^Aj. (1.7) 

i=i 

For the special choice of a staircase partition, 

A = (5„ = (n,n-l,n-2, 1); i{6n) = n, (1.8) 

the identities (1.5),(1.6) and (1.4) lead to the representation of Qn{x) in terms of Schur polynomials. 


Qn{x) = 




(2fc)! 

2f^k\ 


X, 0, —0,0,... I , X G C, nGZ>i. 


It is well known that equation (1.1) admits higher order generalizations and itself forms the first member of 
a full hierarchy. To be more precise, let Cn denote the following quantities expressed in terms of the Lenard 
recursion operator, 

^£Ar+i[u] = + 2u:^ /IatM, N g Z>o; £oM = (1-9) 

and with the integration constant determined uniquely by the requirement £„[0] =0, n > 1. The recursion 
gives, for instance, 

Ci[u]=u, C2[u] = Uxx + , C-i[u] = Uxxxx + ^{Uxf + IQuUxx + ■ 

The A^-th member of the Painleve II hierarchy is subsequently defined as the ordinary differential equation 

Cm\ux — = XU + aM, X G C, Oat G C; u = u{x]aN,N). (1-10) 

Hence, the first member = I is Painleve II (1.1) itself, and more generally, the A^-th member is an ordinary 
differential equation of order 2N. Besides (1.10), we shall also consider a case which involves additional 
complex parameters ts, , t 2 N-i- With u = u(x; aN,t, N) for x, ciAf G C and t = (fa,... , t 2 N-i) G 



/^ci\ 2Vi /^ci\ 

vdx ^ ^ X! l)^ 2 fe+l + 2uJ Ck [ux - +XU + UN- 


( 1 . 11 ) 


For (1.10) and (1.11), it is known [16, 14] that rational solutions exist if and only if = n G Z. Moreover, 
Clarkson and Mansfield in [8] introduced generalizations of the Vorob’ev-Yablonski polynomials for N = 2,3 
which allow to compute the rational solutions of (1.10) once more in terms of logarithmic derivatives, 


u{x; n, N) = — In • 
dx 


Q^u-ljx) 

qL^'(x) 


n G Z>i; it(x; 0, N) = 0, u(x; —n, N) = —u(x; n, N), n G Z>i. 


This approach has been extended to (1.11) for general N G Z>i by Demina and Kudryashov [13, 14] who 
found in particular the analogues of (1.3) for, what we shall call generalized Vorob’ev-Yablonski polynomials 
(x;t). 


dnMx;t)Q\^'‘\{x;t) = (Q[fl(x;f))^<{ x - 2£ 


i[w] 


N 


2^ In Ql^\x;t) 


( 1 . 12 ) 


N-1 


+2 (2k + l)f2fe+i£fc 


k=l 


2^ In Q^n\x-,t) 


>, n G Z>i 


with Q[)^^(x;f) = 1 and Qi^'(x;f) = x. For fixed t = ... ,t 2 N-i) G ^ and n,N G Z>i these 

special polynomials are then used in the construction of the unique rational solutions of (1.11), 


u{x; n, t, N) = — In ■ 
dx 


Q^n-lix;t) 

Qlf^'(x;f) 


u{x; 0, t, N) = 0, u(x; —n, t, N) = — m(x; n, t, N). 
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Figure 1. The roots of the rescaled higher Vorob’ev-Yablonski polynomials Qgg, 

Q 77 , Q 78 > Sts (from left to right and top to bottom). See (1.18), (1.19) for their definition. 
The symmetry of the pattern is easily explained from the definition of the polynomials. The 
locations of the outer vertices of the star shaped regions are given in (1.23). The various 
lines that appear in the Figures are not straight lines but real analytic arcs defined by the 
implicit equation (4.20). It is quite evident that for N >2 there are further subdivisions of 
the star-shaped region into subregions. 


1.1. Determinantal identities. It is mentioned in [13], but not proven, that also Q[f^^(a;;t) can be ex¬ 
pressed as a Schur polynomial. In our first Theorem below we shall close this small gap. 


Theorem 1.2. Let Sn denote the staircase partition (1.8) of length n S Z>i. For any 

t = (ta, ts,..., t2N-i) € 

the generalized Vorob’ev-Yablonski polynomial Q\^\x;t),x G C defined in (L12) equals 

= n ^ 0 . 22^3, 0, 2 %, . . . , 2 ^^t 2 N+l, 0 , 0, 0, . . .) , t2N+l = " ' (^l^) 

Besides the Jacobi-Trudi type identity (1.4), Vorob’ev-Yablonski polynomials can also be expressed as 
Hankel determinants, in fact in [5] the following Hankel determinant representation for the squared polyno¬ 
mial Qn{x) was obtained. 








■(2fc)!]^ 

“fcT 


det [ni+j- 2 (x)]'^^l^, 


X G C 


(1.14) 
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with {/4fc(a;)}fcgz>o defined by the generating function 

31 oo 

exp xw -— = fj,j{x)w^. 

O 

L J i=0 

In our second Theorem we present the analogue of (1-14) for the generalized Vorob’ev-Yablonski polynomial 
Ql^\x-,t). 

Theorem 1.3. Lett = {t^,... ,t 2 N-i) G andn G Z>i. For any x G C we have the Hankel determinant 

representation 


where we use the abbreviation 




(2fc)! 

kl 


det 


4+i-2(to) 


1 n+1 




(1.15) 


to — (^ 1 ,0, ta, 0, ts,..., 0 ,t 2 Ar+i, 0, 0,0,...); ti — x, t 2 Ar+i —— 

and the coefficients defined by the generating function 


1 


27V + 1 


exp 


7=1 




T- =0, j>2N+ 1. 


(1.16) 


k=0 


Remark 1.4. In fact, the statement of Theorem 1.3 is the specialization of a more general identity for Schur 
functions (compare Lemma 3.1 below) which in our case reads 

sl^ (ti, 0,2^73,0, 2%, 0,...) = 2"%(„+i)„ (ti, 0, t3,0, is, 0, tr,...). (1.17) 

Here, A = (n + 1)" denotes the rectangular partition with n + 1 rows of length n and the specialization consists 
in simply setting 

ti = X, t2N+i = ~ _|_ I ’ ^ ^ 

1.2. Roots of higher Vorob’ev-Yablonski polynomials. In analogy to [5], we provide a direct appli¬ 
cation of Theorem 1.3. Numerical studies carried out in [19, 8, 14] show that the zeros of generalized 
Vorob’ev-Yablonski polynomials form highly regular and symmetric patterns as can be clearly seen in Figure 
1. These patterns in case of the Painleve II equation itself have been first analyzed in [6, 7]. However, the 
approach outlined in [5] starts directly from (1.14) and not from a Lax pair associated with (1.1). To be more 
precise, the identity (1.15) allows us to localize the roots of the generalized Vorob’ev-Yablonski polynomials 
as n —>■ oo by analyzing associated pseudo-orthogonal polynomials. Of course in the generalized case these 
patterns depend on the parameters (compare [14]); we shall confine ourselves here to the case of 

higher Vorob’ev-Yablonski polynomials, namely the case 

O = ^5 = ■ • ■ = t2N-l = 0. (1-18) 

More specifically, we are considering the roots of the rescaled higher Vorob’ev-Yablonski polynomials 

fH{fl = |a:eC: Qlf'l = o|. (1.19) 

These sets admit a discrete 'Z 2 N +1 rotational symmetry, which follows immediately from Theorem 1.2 and 
the homogeneity (1.7); 

Q[)^'(a;x) = ui = . (1-20) 

We can provide a partial analytic description for the boundary of the polygons Pn seen in Figure 1 which 
asymptotically contain the sets (1-19) as n —>■ 00 . More precisely we have first the following Theorem. 


Theorem 1.5. There exists a compact region Pn in the complex x-plane such that for any ^ > 0 the root 
sets are contained in a 5-neighborhood JVs{Pn) of Pn as n —>■ 00 . 


The description of the regions Pn is provided in part by Theorem 1.7 below. First we require 
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Definition 1.6. Given N S Z>i, let a = a(x; N), x G C denote the unique solution of the algebraie equation 

{2af^+^ - x{2af^ + (-1)^ (T) = ° 

which is analytic in the domain 


2 N 

X e c\ U [' 




and behaves near x = oo as 


Here, the points x = x^^\ fc = 0,..., 2N are the solutions of 


a = — + O (x , X —>■ oo. 


..27V+1 


= (-1)^(27V+1) 


/27V + /27V 


V 2N 


\ N 


( 1 . 22 ) 


(1.23) 


and form the outer vertiees of the regular star-shaped regions shown in Figure 1. 

Theorem 1.7. The regions Pm are compact, invariant under the rotations of angle 2 W^’ contain the origin 
and their boundary OPm consists of branches of the locus in the complex x-plane described by 


Bn = ja:: S C : 3 ?((/j( 2:; a)) 


AN] 


= 0^ 


(1.24) 


Here z = k = 1,..., 2N are the solutions of the equation 


whereTm,a{z) denotes the Maclaurin polynomial of degree m € Z>o of the function (H-z)“ = l-\-0(z),z —>■ 0. 
Moreover a = a(x) is defined in (1-21) and (1.22), and the function ip is defined by 


(p{z; a) = —2 In 


2: + (z^ + a^) 2 


+ -(z^ + 0 ^)" - 


1 (z2+a2); 


T., 


2iV + l a3z2^+i 


with principal branches for fractional exponents and logarithms. 

The branches of the real-analytic curves specified by Biv of Theorem 1.7 are plotted as the arcs in Figure 
1. Perhaps more important than what Theorem 1.7 above says, is what it does not say. In fact of all the 
branches of curves defined by (1.24) we are not able to effectively discern which ones actually form the 
boundary of Pm- In particular we cannot conclude in general that the points (1.23) belong to OPm- 

Remark 1.8. A local analysis (which we do not propose here but is essentially identical to [5]j shows that 
the angles between consecutive ares emanating from the points x^^^ (1.23) is 

1.3. The roots inside Pm- Inspection of Figure 1 clearly shows that the pattern of roots within Pm 
is subdivided in subregions. This can be easily qualitatively understood in terms of the steepest descent 
analysis; the so-called g-function of the problem (see Section 4 below) is an Abelian integral on a Riemann 
surface of genus 0 on the outside of Pm and of genera 2,4,... inside. In fact we can show that x = 0 
belongs to a region where the genus is 2N and thus it is reasonable to deduce that there are nested regions 
of higher and higher genus, until the maximum is reached (2N). These regions are quite evident in Figure 
1. In principle the boundaries between these nested regions could be described as well in terms of Abelian 
integrals, but it is beyond the scope of this paper to attempt any such detailed description. 
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1.4. Outline of paper. We conclude the introduction with a short outline of the upcoming sections. First 
Theorem 1.2 is derived in Section 2 by referring to the KdV and mKdV hierarchies for which we construct 
a rational tau function in terms of Schur polynomials. Subsequently an explicit scaling reduction brings us 
back to the Painleve II hierarchy and Theorem 1.2 follows. After that we turn towards Theorem 1.3, but 
opposed to the proof of (1-14) in [5] which relied on (1.3), identity (1.15) will follow from Schur function 
identities and Theorem 1.2. In the final Section 4 we follow largely the logic outlined in [5]. However we 
choose not to present any details on the nonlinear steepest descent analysis for the underlying orthogonal 
polynomials. Once the correct inequalities for the g-function have been verified the asymptotic analysis 
outside of Pn is almost identical to [5], see Section 4 for further details. 

1.5. Acknowledgments. All authors are grateful to P. Clarkson for useful discussions about this project. 
M.B. is supported in part by the Natural Sciences and Engineering Research Council of Canada. F. B. and 
M. B. are supported in part by the Fonds de recherche Nature et technology du Quebec. T.B. acknowledges 
hospitality of SISSA, Trieste in February 2015. Early stages of the manuscript were carried out while F.B. 
was a Postdoctoral fellow at SISSA. 


2. Short reminder about the (m)KdV and Painleve II hierarchies 

The goal of this section is to remind the reader very briefly of the construction of the Painleve II hierarchy 
as a scaling reduction of the modified Korteweg-de Vries (mKdV) hierarchy, cf. [8]. In doing so we will en 
route derive Theorem 1.2. 


2.1. The KdV hierarchy. The KdV hierarchy involves the Lenard recursion operator 
d { d \ 1 

—Cn+i[u\=i-^+'iu{x)—+2u^{x)jCn[u\, Tq M = -, [0] = 0 

and its equations are written as the partial differential equations 
du d 

— -= — £„+i[m], nGZ>o; u = u{to), = (ti, 0, fa, 0,^5,...). 

dt2n+l dx 


( 2 . 1 ) 


( 2 . 2 ) 


It is customary, and we will adhere to the custom, to denote the variable ti by x since Li[u] = u and hence 
the first member of the hierarchy above reads simply dt^^u = Ux- In general, the equations of the hierarchy 
should be viewed as an infinite set of compatible evolution equations for a single function u = u(x). A 
solution of the hierarchy is then a function u{x; fa, ts,...). 


Definition 2.1. A function = Tj^dv(^o) called a tau function for the KdV hierarchy (2.2) if the 
function 


solves the hierarchy (2.2). 
tau function. 


u{to) = 2^1nrj,dv(to) 


We note that multiplication by an 


, x = ti. 

arbitrary constant (in x) of 


gives another 


2.2. Rational solutions to KdV and staircase Schur polynomials. The solutions of the KdV equation 
rational in x for all values oit = t^ (and for all higher times ts, ty,...) and vanishing at x = oo were completely 
characterized in [2]; they all belong to the countable union of orbits flowing out of initial data of the form 

ti(ti l ) 

u„(x, 0,0,...) = —- , n G Z>0 . (2.3) 

The corresponding tau functions r„(to) were obtained explicitly in [1] in terms of Wronskians of certain poly¬ 
nomials in tp. Up to normalization and re-parametrization these Wronskians coincide with Schur polynomials 
associated to staircase partitions evaluated at the odd times, namely 

^2 

Tnito) = Ss„iti,0,2‘^t3,0,2%,...) , u„(to) = 2—lnr„(to) , (2.4) 

where d„ denotes the staircase partition (1.8) of length n G Z>i. Moreover, it can be shown (cf. [21]) that 
these are the only Schur polynomials that give KdV tau functions when all even times are set to zero. 
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Remark 2.2. The particular rescaling >—>■ is used in (2.4) in order to correct the normalization 

so that the coefficients in (2.2) are as indicated. 


2.3. The mKdV hierarchy. The modified KdV (mKdV) hierarchy is defined in terms of a new dependent 
variable v = v(to) which is related to u via the Miura transformation 

u = TVx-v‘^, (2.5) 

where the choice of signs is arbitrary. More is true: if v satisfies u = —Vx — then the new function 
u = Vx — v'^ is a different solution of the KdV hierarchy (and vice versa); this is an example of a Backlund 
transformation. Inserting (2.5) into (2.2) yields a new set of evolution equations 


d 


dt 


2n+l 


V ^ ^ Vx f ^ ^ '^{^Vx H“ 2(dz'Ua;a; 2vVx) J ^ Vx V 


This can be rewritten as follows 

d „ 
ax 


dv 


dt 


2n+l 


or equivalently 


d „ 

- 2n 
ox 


d_ 
dx 

dv d 
dt2n+i dx 


^r.[TVx-v^] 1=0. 


( 2 . 6 ) 


h] 


We now notice that the two expressions 
d 




dx 


da: ( ^ ± 2u ) ^ T 2u 
dx ) dx \dx 


(2.7) 


define the same differential polynomial in v since the right hand side is clearly invariant under the map 
x i-A —X. Thus we can simply write 

£3WH = O1-)H=0„H, 

omitting the reference to the choice of sign. We now want to conclude that the expression 0ri[n] vanishes 
identically; the two equations in (2.8) below are simply stating that F{x) = 0n[f] is a joint solution of 
the two ordinary differential equations (Fdx + 2v)F{x) = 0. Thus 0„[i;] should be the in the null-space of 
both equations Fdx + 2t>; as long as v is not identically zero (which is an un-interesting situation), the only 
function in both null-spaces is the null function and hence 0„[u] = 0. Thus we have concluded that if u 
is a solution of the KdV hierarchy (2.2) and v is related to u by (2.5), then v must solve the hierarchy of 
equations indicated below and named mKdV hierarchy, 

du d f d X 

dff —^ ~ 71 G Z>o; v = v(to). (2.8) 

The choice of signs is irrelevant, since the right hand side (as noted above) yields the same differential 
polynomial in v. 


2.4. Schur functions and Painleve II hierarchy. Let us now return to our special situation for which 
we fix 

ti=X, t2N+l = ~ 2_/y _|_ ^ = (^3, ^5, ■ ■ ■ , i2N-l) G t2j + l = 0, j>N. 

Proposition 2.3. For n,N G Z>i define the two functions 

gn{x;t) = In S5^(a:,0,2^t3,0,2'‘t5,...,2^^t2iv+i, 0,0,0,...), Wffxff) = gn+iix;t) - gn{x]t) (2.9) 
with some fixed branch for the logarithm. We then have the Miura relation 

2dlgn{x-,t) = -dlWn{x;t) - {dxWn{x;t))^. 


( 2 . 10 ) 
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A proof of (2.10) can be found in Appendix A. In view of Proposition 2.3 we note that the two functions 

u(to) =2dlgn{x\t), v{to) = dxWn{,X\f) ( 2 . 11 ) 

satisfy precisely the Miura relation (2.5) with the choice of the minus sign, namely u = —v' — v^. Since ss^ 
gives a tau function for the KdV hierarchy it follows that v satisfies the hierarchy (2.8) for n = 0,... ,1V. 
Summarizing 


Proposition 2.4. The function 

satisfies the mKdV hierarchy in the form 

dw 


j{x;t) = -dxWn{x;t) 


dt 


2n+l 


= ^ + 2ui^ Cn [wx-w^], 0 < n < N. 


( 2 . 12 ) 


Recalling the homogeneity property (1.7) we see that w{x]tf) obeys a simple scaling invariance which will 
allow us to reduce the partial differential equations (2.12) to an ordinary differential equation; we carry out 
a scaling reduction: 

(i) View w = —i;(to) as a function in the variables ti = x,t = {t^^t^,... ,t 2 W-i) G and t 2 N+i- 

(ii) By homogeneity (1.7), it follows that w = w(ti,t,t 2 N+i) is a function of the form 


W={-{2N + l)t2N+l) ^~ + ^P(ri,r3,...,T2Ar_i), 
and V depends on the “new” variables 

T2k+i = - — -wr, k = 0,...,N-l. 

{—{2N + l)t2Ar+i) 

(iii) Substituting (2.13) into the left hand side of (2.12) with n = N, we find 

IV-l 


(2.13) 


dw 


{2N + l)t2N+i-^ - = ~ (2-^ + l)i2Ar+i) 


V + ^ (2j + 1)T, 
1=0 


dV 


2j + l gj. 


21 + 1 


(iv) Next we evaluate (2.14),(2.13) at t 2 N+i = — 2 n+i ^’^ol compare the result to (2.12), 


d^fd 

dx 


f 3 \ 

f -g —I" 2V j Ln [14 — = V + (2j + l)t2j+i 

V ^ / j=o 


dV 


dt 


2j + l 


(2.14) 


(2.15) 


(v) Since ti = x and V + = ^{xV), (2.15) can be rewritten with the help of (2.12), 


dx 


{ic + -V^]-xV-^ (2j + 1)4 ,+1 + 2v) C, [V, -V^]\= 0. (2.16) 


Equation (2.16) is an ordinary differential equation for the function V = w{x',t) in which t G ^ appear 
as parameters. Since 

TL I 1. 

w{x-,t) = da;{gn{x;t) - gn+i{x;t)) = - - -h O (x~^) , x ^ 00 , 

integration in (2.16) yields (1-11) with ajy = n + 1. Recall [13, 14] that aj^ G Z is necessary to have a 
rational solution to (1.11) and for all integer values oi a n there exists a unique rational solution which 
can be obtained from the trivial solution for a^v = 0 by Backlund transformations. Therefore we have the 
following 

Theorem 2.5. For n,N G Z>i the unique rational solution of the Painleve II hierarchy (1.11) is 

u{x;n + l,t,N) = ^ In (x, 0, 22t3, 0, 2^4,..., 22^4Ar+i, 0, 0, 0,...), t 2 iv+i =(2-17) 


dx 


2V+ 1 
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and we have the identity 


Q\^\x;t) = ^^^ss^{x,0,2'^t3,0,2%,... ,2'^^t2N+uO,0,0,...), x € C. 


k=l 


Proof. It is easy to see that the LHS of (2.17) is a rational solution to (1-11) by the scaling reduction 
(2.13)-(2.16). By the uniqueness of the rational solutions of the Painleve II hierarchy we have 

t) = c„,Ar(f)s 5 „ {x, 0, 2 ^t 3 , 0, 2%, 2 ^^t 2 Ar+i, 0, 0,0,...), 

with an a;-independent factor Cn,Nit). However, to leading order, 

ss^ (x, 0, 2%, 0,2%,..., 2^^t2N+i,0, 0,0,...) ~ S5„ (x, 0,0,0,...) = 




where h{X) denotes the product of the hook-lengths of A (cf. [20]). Since 


O ^ I 

l<5n| = ^(n+l), /,(d„) = Jl — 


and Q.\^\x;t) is a monic polynomial of degree ^{n+ 1), the claim follows. 


□ 


3. Proof of Theorem 1.3 

We will appeal to certain identities satisfied by symmetric functions which can be found, for instance, in 
[20]. First let us start with the following lemma. 

Lemma 3.1. The symmetric polynomial identity 

(2°ti, 0, 2h3, 0,24^5,0,...) = 2-"s(„+i). (2ti, 0, 2 H 3 . 0,2%,...) = 2"%(„+i)„ (ti, 0, h, 0, h, 0, tr, • ■ •) 

holds, where X = {n + 1)" stands for the rectangular partition with n -|- 1 rows of length n and 6n is the 
staircase partition (1.8). 

Proof The Schur polynomial 5 , 5 ^ can be written in terms of the projective Schur polynomial Ps^ labeled by 
the same partition 

S5„(tl,0,t3,0,t5,0, ...) = Ps^{ti,t3,t5,...). (3.1) 

For a proof of this identity see [20], § 3.8, example 3, page 259, and also [15], Lemma V.4. Second, for a 
strict partition A, i.e. Ai > A 2 > ... > A^(;i,), we have [25], Theorem 4, 


(|, I’ I’ ■ • •) = ®A(il, 0, t3, 0, is,...), 


(3.2) 


with A denoting the double of the partition A which is defined via its Frobenius characteristics, 

A = (Ai, A 2 ,..., A^(a) I Ai — 1 , A 2 — 1 ,..., A^(a) — 1 ). 

Combining (3.1) and (3.2), 

( 2 °ti, 0 , 2 h 3 , 0 , 24 ^ 5 , 0 ,...) = Pi { 2 ^, 2 %, 2 %,...) = 2 -”s(„+i)„ ( 2 ti, 0 , 2 H 3 , 0 , 2 %,...) 

= 2" S(„+l)n(tl, 0, t 3 , 0, ts, ...) 

where we have used homogeneity (1.7) in the last step. This concludes the proof. 

We are now ready to derive Theorem 1.3 by referring to (1.13) and Lemma 3.1. 

Proof of Theorem 1.3. Let t= (^ 3 , ts,..., t 2 Ar-i) G and 

1 


□ 


to = (ti, 0, ^3, 0, is, . . . ,f2Af-l, 0, <2Ar+l, 0, 0, 0, . . .), ti=X, t2N+l= — 


2N + I' 
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This gives us 


n 

fc=i 

1 

2 " 

1 


(2fc)! 

2>^k\ 


( 2 %, 0 , 2%, 0 , 2%,2^^t2N+uO, 0 , 0 ,...) 


n 


n 


M! 

ki 

im 

k\ 


n r \ n 2 


S(n+l)"(il,0,t3,0,t5, . . . , t2Ar+l, 0, 0, 0, . . .) 

n r Ml 2 


<L )=4 n T 


'^(n)^+i i 




^ r /r^ 7 \ n 2 


A;=l 


fc=l 

[N] 


(2fc)! 


where we used that for the transposed partition A', 

SA'(ti, 0, ta, 0, ts,...) = (—1 )''''sa(— ti, 0, —ts, —ts,...) ^ = ^ SA(ti, 0, ta, 0, ^5,...), 

and that the Schur polynomials of rectangular partitions are Hankel determinants. 

Corollary 3.2. Let to = (ti, 0, ta, 0, ts,...) and {hk{to)}kez>o (^■^)- Introducing the notation 

^n,I!(to) — det (to)] ; n, £ € ^> 0 , 

we have the Hankel determinant identity 


Ara+l,o(to) — ( —l)"A„_2(to)- 


(3.3) 

□ 


(3.4) 


Proof. Note that 

'S(ra+i)" (to) = det [/in+i—j+fe(to)]^-^^_j^ = (~ 1 ) det [hj_|_fc(to)]^-^^_j^, 

S(„)..+i(to) = det [/i„_j+fe(to)]"j];]^^ = (-!)"■ det [/rj+fe_ 2 (to)]”^;]^p 

and since |(n)"+^| = |(n + 1)"| = n{n + 1) = Omod 2, the stated identity follows from (3.3). □ 

Remark 3.3. Identity (3.4) in Corollary 3.2 does not hold if any of the even-index times is nonzero. 


4. Characterization of the set 

The logic we are following here is identical to [5] . The square of the polynomials (x) is proportional 
to a Hankel determinant 

An{x;N) = det [Aij+i_ 2 (to)]”tii 
of the moments fJj^\to) (1.16), which can alternatively be written as 

Hk i'^o) = to = (x, 0 , 0 , ..., 0 ,t 2 iv+i, 0 , 0 , 0 ,...), t 2 jv+i = ~ 2 jv + l 

where S' C C denotes the unit circle traversed in counterclockwise direction. It is then a well-known fact 
that An{x; N) = 0 if and only if the Riemann-Hilbert problem 4.1 has no solution, or equivalently, if and 
only if the n-th monic orthogonal polynomial for the weight 

g-i9(z;x,Af) 2 X 

dMo(^;x,iV) = ^-—, + ^es 

does not exist. In view of the scaling x 1 —>■ n^^+ix in (1.19) we also perform a scaling z 1 —>■ z so that 

we arrive at the following Riemann-Hilbert problem with a varying exponential weight. 

Riemann-Hilbert Problem 4.1. Suppose 7 C C is a smooth Jordan curve whieh eneircles the origin 
counterclockwise. Let T = T[z',x,n,N) denote the 2x2 matrix-valued piecewise analytie function which is 
uniquely characterized by the following three properties. 

(1) r(z) is analytic for z G C \7 
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(2) Given the orientation of^, the limiting values r±(z) from the (+) and (—) side of the contour exist 
and are related via the jump condition 

Q—n'&{z;x,N) 


r+(z) = r_(z) 


1 w{z]x,N) 
0 1 


z e 7 ; w{z] X, N) = 


2TTiz 


(3) The function r(z) is normalized as z ^ oo, 

Ti{x;n,N) 


m= 1+ 


^ + 0 (z-2) 


^3 = 


1 0 
0 -1 


Then we have, compare [5], 

Proposition 4.2. The zeros of the scaled Vorob’ev-Yablonski polynomials (n^w+i x) coincide with the 
values of X for which the problem j.l is not solvable. 


In principle an asymptotic analysis of the Problem 4.1 as n —>■ oo is possible using the Deift-Zhou steepest 
descent analysis [9, 11, 12], and the zeros will be located asymptotically in terms of appropriate Theta 
functions as in [5]. However here we simply want to prove the absence of zeros outside of a certain compact 
region Pjv and, en route., give a partial characterization of the boundary dP^. For a more comprehensive 
analysis which is only marginally different from the present situation we refer to [5]; here we shall just remind 
the reader that the method requires the construction of an appropriate function, called customarily ’’the 
0 -function”. 

In case of the problem 4.1 the g-function is a priori expressible in the form 

Qiz) = ^'diz;x,N) + J (P 4 jv+ 2 (w))" ^ 2 N +2 + i = iix;N)GC, zq = zo{x; N) G C (4.2) 

where, in general, P 4 N+ 2 iz) is an appropriate polynomial of the indicated degree. The ansatz (4.2) is 
explained in the paragraph ’’Construction of the g-function” of [5] and the discussion there can be applied 
almost verbatim here. From (4.2) we see that the g-function is an Abelian integral on the Riemann surface 
of the square root of P 4 .n+ 2 {z); depending on the number of odd roots, this surface has a genus that can 
range from a minimum of 0 (if there are only two simple roots in P 4 N+ 2 iz)) to a maximum of 2N (if all the 
roots are simple). Subsequently the Deift-Zhou analysis shows that 

If X is such that the genus of the above Riemann surface is zero and the g-function satisfies the 
appropriate inequalities (recalled below), then the Riemann-Hilbert problem 4.1 is eventually 
solvable for sufficiently large n. 

Therefore our strategy is as follows; we postulate a genus zero Ansatz for the g-function in (4.3); the algebraic 
requirements are easily verified, but the required inequalities are not always verified. We shall then find where 
the inequalities fail, and hence where the roots are asymptotically confined. 

For the concrete construction of the g-function in the genus zero region we follow the logic outlined in [5] . 
We seek a function y = y{z), z G <C\B of the form 

= PeC[z] deg(P) = 2iV, (4.3) 

where (z^ -I- a^)5 is defined and analytic off the oriented branch cut B = B{x, N) which connects the points 
z = ±ia. The precise location of P C C shall be discussed in Section 4.2 below, for now we require that y 
satisfies the two conditions 

y{z) = ^'&ziz;x,N)+ 0{1), z^O- y{z) = ^ + O {z~'^) , z ^ oo. (4.4) 

Using simple algebra, we directly obtain 


Proposition 4.3. The conditions (4-4) imply that a and x are related via 

„ (-1)^ 


Cn 

Pm' 


22N 


X = 2a-\- 


Cn = 


2N 

N 


(4.5) 
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and the polynomial P{z; a) is uniquely determined as 




(4.6) 


where Tm,a(,0 is the Maclaurin polynomial of degree m S Z>o of the function (1 + C)“ = 1 + 0{f), C —)• 0. 

Proof. Observe that the condition y{z) ~ z —>■ oo implies that P{z; a) is monic and from the behavior at 
z = 0 we find that 

= + (4.7) 

z V2 + a 

Writing P{z-,a) = z'^^ + Q{z) with a polynomial Q{z) of degree at most 2N —1 and reading (4.7) at O (z'^^), 
we get 

= = Q{z) + O {z^^) , z —t 0, 


2^z 


and thus 


(;;2 


which gives (4.6). In order to deduce (4.5), we recall that P{z;a) is monic, i.e. we must have 


— ( ^ 
j2A7 1 ^ 


= 2a, 


that is (4.5). 


□ 


Proposition 4.4. The branch points of the map a = a{x), defined implicitly by (4-5), coincide with the 
values of a{x) for which two zeros of P(z) overlap with the branch points z = ±ia. 


Proof. We have to evaluate the condition P{±ia) = 0; using (4.6) this amounts to 


N2N 


0 = P(±ia;a) = (-l)'''a 


--V 

2a ^ 


N-l / 1 

2 


fc=0 




\N 


N ( - 


a \ N 


N2N ^ 


= (_l)iVa^iV _ 


a22^ \ N )' 

Thus the condition determining the coincidence of a zero of P{z) with 2 = ±ia is 

2A+1 _ (-1)^^ /^2iV' 


22N y JV 

On the other hand the map (4.5) has a branch point where x'{a) = 0, which gives exactly (4.8). 


(4.8) 

□ 


4.1. The complex effective potential and the inequalities. For further steps it will prove useful to 
define the effective potential, 


ip{z\a) = 2 / y{w)dw, z S C\S 

J laix') 


%{x) 

which in the given situation (4.6) can be evaluated explicitly. 


ip{z] a) = —2 In 


z + (z^ + a^)2 


1 


10 




-T.; 


(z^ + a^)l 

q3^2A7+1 


(4.9) 


(4.10) 


and all branches in (4.10) are principal ones such that (z^ + 0^)2 ~ z as z —>■ 00. 
Lemma 4.5. Given ( 4 . 10), we have for N G Z>i, 

^ (z^+a^)5 


(4.11) 
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Proof. The jump of (p{z; a) equals 47ri on a contour that extends to infinity. Hence dxf{z', a) has no jump on 
a contour which extends to infinity. Along B we have {dx(p{z; a))+ = —{dx^p{z; a))-. Since (p{z;a) vanishes 
at the branch point z = ia (and is constant ±27ri at z = —ia on the two sides) we deduce that dxf{z;a) 
must be zero at ±ia. Also (compare Section 4.2 below), 


z^O 

dx^ ’ ’ | 0 ( 1 ), 


(4.12) 


Thus the ratio of the proposed expression for d^^iz; a) is bounded at z = ±ia, analytic across the cut and 
bounded at z = oo with limit 1 at z = 0. The Proposition now follows from Liouville’s theorem. □ 


The potential (4.9) is related to the g-function (4.2) by 

Siz) = ^{t}{z;x,N) - ip{z;a) + i), z G C\B (4-13) 

where the constant £ {modified Robin constant) is defined by the requirement that g(z) = ln(z) + 0(z~^) as 
|z| —>■ oo. The relevant inequalities for g(z) are more conveniently expressed directly as inequalities for the 
effective potential. In terms of the latter, the following properties of the effective potential are equivalent to 
the existence of the g-function and characterize <p(z; a) (the proof of these statements is simple if not already 
obvious) 

(1) Near z = 0 the effective potential has the behavior 

<p(z;a) =-'d(z;x) + 0(1) p(z) = ^i9;,(z; x) + 0(1) , (4.14) 

while near z = oo it behaves as 

(/j(z) =-21nz-tC>(l). (4.15) 

(2) Analytic continuation of <p(z; a) in the domain z G C\B yields the same function up to addition of 
imaginary constants; in particular, the analytic continuation of (p(z; a) around a large circle yields 
(p(z-, a) + 27ri; 

(3) For each component Bj of B we have that, 

93 +(z; a) + P-(z; a) = —2iaj, z G Bj, aj G K 

(4) The effective potential 

$(z; x) = 3?(i^(z; a)), z G C\B (4-16) 

with a = a(x) as in (1.21) and (1.22) is a harmonic function in z G C\B. Moreover <i)(z; •)|g = 0. 

(5) Inequality 1. The sign of 4>(z) on the left and right of B is negative. 

( 6 ) Inequality 2. We can continuously deform the contour of integration 7 to a simple Jordan curve 
(still denoted by 7 ) such that H C 7 and such that <i>(z)|^^g > 0. 

Note that (p(z; a) and g(z) are both related to the antiderivative of the differential 

i (\z 

y(z)dz = (P4 n+2(z))^ 

which is defined on a Riemann surface X = {(ic,z) : w'^ = P 4 n+ 2 (z)} of genus between 0 and 2N. Since 
<i>(z) = 3?((^(z)) vanishes along B, it also follows that H is a subset of its zero level set; therefore, B consists 
of an union of arcs defined locally by the differential equation 3?(?/(z)dz) = 0. 

4.2. Location of branch cut. The following Proposition appeared in [5] but applies also to the present 
situation (4.10). 

Proposition 4.6. The effective potential $(z;a) = 5 ft((/ 3 (z;a)) has the following properties: 

(1) The function $(z;a) is defined modulo a sign depending on the determination of (z'^ + a'^)^. 

(2) The zero-level set Z = {z G C : 4’(z; a) = 0} is well defined independent of the determination of the 
square root in (1) and invariant under the reflection z 1 —> —z. 

(3) For |a| sufficiently large there are two smooth branches of the zero-level set Z which connect z = ±io 
and which are symmetric under z 1 —>■ — z. 
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Proof. Statements (1) and (2) follow just as in 

/ ^ , oi f z + {z^ + 

(p{za;aj —>■ —2 In — ^-- 


Proposition 3.5), for (3) we note that, as a —>■ oo 


'(^2 + 1)2 = _2 




* - 0(4 


(4.17) 


and the limit is uniform on compact subsets of the Riemann sphere not containing z = 0. The remaining 
logic is now as in [5]. □ 


Suppose that a > 0 is sufficiently large and thus Proposition 4.6 applies. We claim that B = B{x, N) is 
the branch in point (3) above which intersects the positive half ray K+ (by deformation this fixes the branch 
cut for all X ^ Pn)- In order to see this, recall from (4.14), as z —)■ 0, 

viz) -^d^(z;a:,TV) ^ 2^27V+2 ^ y{z) ^+oo, zfQ. 

But this requires in (4.3) that (z^ + a ^)5 ~ —a as z —>■ 0. Simultaneously (4.15) requires y{z) ~ \ near 
z = oo and hence (z^ + a^)5 ~ z as z —>■ + 00 . Hence, the determination of the square root in (4.3) has to 
change on the positive half ray, i.e. B = B{x,N) is as claimed. 

4.3. The inequalities of the g-function and the region P^. Since the quadratic differentials 

dz^ dz^ 

ry= (d<^)2 = 4 (z 2 +a^)P^(z;a)^ 4 j^, and in general y = {dip)^ = AP4N+2{z)^^j;^, (4.18) 

are of the type studied by Jenkins and Spencer [17], that is, of the form i?(z)dz^ with i?(z) a rational 
function, we can follow some of the reasoning which was already explained in [5]. 

Preliminaries. Define the set to consist of the union of the second order poles and all critical 
trajectories, i.e., all solutions of 3fi(di^(z; o)) = lft( 2 ?/(z)dz) = 0 that issue from each of the zeros and simple 
poles of i?(z); the latter are absent in our case. The zeros are at ±ia and at the N pairs ±Zj which are the 
roots of the even polynomial P{z;a). Also [22], there are 2fc + 1 branches of issuing from each of the 
points of order k of R{z), k = —1,0,1,... (the case fc = —1 corresponds to simple poles, and all others to 
zeroes). We are interested in the connected components of 

3 

and a simple argument in analytic function theory (see [17]) shows that each simply connected component 
Kj is conformally mapped by (p{z;a) into a half-plane or a vertical strip a < <l>(z;a;) < fi] each doubly 
connected component K is mapped to an annulus (or a punctured disk) {r_ < |?w| < r+} hy w = 
where p = and 7 is a closed simple contour separating the two boundary components of K. It is 

also shown in [17] that there are no other possibilities for the topology of the connected components Kj. 
Moreover, there is a one-to-one correspondence between annular domains (including the degenerate case of 
a punctured disk) and free homotopy classes of simple closed contours 7 for which d(p 7 ^ 0 . In our case 
there is only one such class corresponding to a loop encircling the origin, and hence only one annular domain 
which we denote by K^o (which is actually a punctured disk). 

By construction, (p(z;a) ~ —21nz, z —> cx) which shows that z = 00 is at the center of a conformal 
punctured disk via the conformal map w = Moreover the level sets Cr = {z : $(z;x) = — 21nr} are 

foliating a region around z = cx) in topological circles if r is sufficiently large. Thus none of the hyperelliptic 
trajectories issuing from ±ia, {±Zj} can “escape” to infinity; they either connect to z = 0 or amongst each 
other. Suppose rp is the infimum of the r > 0 for which Cr is smooth; this means that Crg contains at least 
one zero of dz(p (by symmetry, it contains then two zeros in our situation). The annular (punctured disk) 
domain is then (see Figure 2) 

= U (4.19) 

r>ro 


We denote also Dq = C \ AToo, which is a simply connected, symmetric region containing the origin. 
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Necessary and sufficient condition for the correct 
inequalities in genus zero. We argue that we need to have 
ro = 0. To put it differently, the “first encounter” of the level 
sets Cr as r decreases must be with the two branch points ±ia 
rather than any of the zeros {zj}. We shall then verify that this 
occurs for x > 0 large enough. 

Sufficiency. Suppose now that tq = 0 and thus ±ia G K^o and 
±Zj G Int(£)o). Then the simple, closed loop is separated 
into two components by ±ia and each of them is an hyperelliptic 
trajectory. We know that there must be three trajectories from 
each ±ia and two of them are already accounted for and form 
the boundary of Dq (see Fig. 2); thus the third trajectory is 
entirely contained in Dq, which is compact. 

Now let us turn our attention to Dq] the points ±Zj G Dq for 
j = 1 ,... ,iV. In Do each branch of y{z) (4.3) is single valued 
(the branch points are on the boundary of Dq). Only one of 
the two branches of y{z) has the behavior ^t}z(z]x)] integrating 
this branch from ia coincides with ip(z]x) in Dq. The value of 
the sign of <i> in the interior of Dq close to the boundary Dq 
determines which of the two parts of i9Do\{±ia} is the branch 
cut B: this is the part which has $ > 0 on both sides (ie. in Dq 
and Koo). Thus $ is continuous but not harmonic on B, while 
on dDQ\B it is continuous and harmonic. We still need to show 
that there is a path connecting ±ia and which lies within the 
region $ < 0 . 

This follows from the topological description of the possible regions Kj discussed in the paragraph ’’Pre¬ 
liminaries” . Indeed let Ki be the region containing the arc dDQ\B where y}{z] a) is conformally one-to-one. 
From the discussion of signs thus far, this is either a half-plane w = <i> < 0 or a strip — e < $ < 0 (the only 
annular domain is iFoo)- The two points ±ia are mapped on the imaginary axis Jftrc = $ = 0; thus there is a 
path connecting (p(ia) to <p(—ia) in the w-plane which lies in the left half plane. The pre-image of this path 
in the z plane connects thus ±ia and $ restricted to the interior points of this path is strictly negative. 

Necessity. If ro > 0 then ±ia G Dq. The trajectories issuing from ±ia all belong to the zero level set of 
4). None of them can connect to any of the zeros {±Zj}, and thus they either connect to each other or to the 
origin. Since the sign of $ changes 2N + 1 times around z = 0, they all must go to the origin and thus there 
is no possibility of deforming the contour of integration so that it contains the branch cut B and avoids the 
origin. 

Sufficient condition for the correct inequalities in highest genus. We work with the same general 
setup as in the previous case. Now the quadratic differential is of the form on the right in (4.18). Suppose 
that P 4 N+ 2 {z) there has all simple roots (the roots come clearly in pairs of opposite signs). We 

claim that a sufficient condition for the fulfillment of the inequalities is that all branch points lie on 
ODq = dKoo- In this case dKoo is broken into AN + 2 arcs (see for example Figure 3). There is only one 
branch of y{z) that behaves as y{z) ^ |'d 2 (z;x) near z = 0; the integral of this branch with base point 
is single-valued in Dq = C\iFoo because the region contains no branch points and the residue of y{z) at 
z = 0 vanishes; this integral then defines ip{z]a) (and $) within Dq. The level curves of 4>(z;x) that issue 
from Qj and do not connect to other branch points must connect to the origin because <I>(z;x) changes sign 
exactly AN + 2 times when going around the origin. The regions where f{z] a) is now one-to-one within Dq 
are AN + 2 half-planes because their boundary has only one connected component. Necessarily in 2N -|-1 of 
them <i)(z;x) < 0 and 2N -|- 1 of them <i)(z;x) > 0. The arcs of bounding the three regions 

where $(z; x) < 0 are the cuts and the other are simply zero level sets separating regions where $(z; x) has 
opposite signs. The possibility of connecting two branch points that are connected by an arc of these level 
sets follows exactly by the same argument used in the previous paragraph. 



Illustration of the conformal punctured 
disk Kao, foliated by the trajectories Cr. 
The complement, Dq, contains the other 
critical trajectories. In this example N = 
2 so there are 4 saddle points inside Dq, 
visible here where the critical trajectories 
intersect at right angles. 

























GENERALIZED VOROB’EV-YABLONSKI POLYNOMIALS 


16 


Occurrence of the necessary/sufficient conditions. By point (3) of Proposition 4.6, for a (hence, 
x) large enough there is a smooth branch of the zero level set of 4>(z;x) that connects ±ia; by symmetry, 
there is another one and thus the third branch of the level set that issue from ±ia must go to zero (we have 
seen that there is no branch that extends to infinity). The remaining roots of P{z;a) all tend to zero as 
a — )■ 00 (which is easily seen from the explicit expression (4.6)). Thus they must fall within the region Dq. 
Then the necessary condition in genus 0 is fulfilled. For the case of maximal genus, the occurrence of the 
sufficient condition is contained in Proposition 4.9. 

The discriminant locus. By discriminant locus we refer to the boundary of the locus in the x-plane 
where the inequalities for the genus zero Ansatz fail. This is the boundary of a region Pjv; from the discussion 
above it follows that the inequalities are preserved under a deformation in x, until failure occurs exactly 
when one of the zeros of P(z; a) falls on the branch cut B connecting ±ia, and hence, by symmetry, one also 
intersecting the opposite arc. We know that this does not happen for |x| sufficiently large and hence the 
discriminant must be a bounded set. 

In order to detect the occurrence of the situation above it is necessary (but not sufficient) that $(zj; x) = 0 
for some j, i.e., one of the saddle points of $ lies on the zero-level set; 

dPN c3N = {a;€C: 3zeC: 3 fi((/ 3 (z; a(x))) = 0, (pzz{z; a{x)) = O} . (4.20) 

The set 3 at is clearly closed, and thus OPn must be compact (since we know already it is bounded). However, 
the set 3 AT is strictly larger than dP/v; indeed it describes the situation where any of the saddle points of 
4)(z;x) intersect any branch of the zero level set; The zero level set contains several branches besides the 
branch cut and hence the set 3Af in (4.20) describes also all these “fake” situations. 

A detailed analysis for arbitrary N seems unwieldy. We shall attempt below only a partial study of the 
case N = 2,3 in Appendix B, where we show that the points (1.23) do indeed belong to the boundary of 
the regions Pn for N = 2,3. However the set 3Af is easily drawn and the results are displayed in Figure 1 
together with the roots of some higher polynomials. The result of this discussion is the following theorem; 

Theorem 4.7. The roots of the polynomials x^ lie all within an arbitrarily small neigh¬ 

borhood of a compact region Pjv as n —>■ oo; the boundary of this region consists entirely of a finite union of 
real-analytic arcs in the x plane satisfying the condition (4-20). 

The condition (4.20) is spelt out in more detail in the statement of Theorem 1.7, which is henceforth 
proved as well. 

Remark 4.8. A careful consideration should allow also to prove that the region is simply connected. It is 
also relatively simple to show that x = 0 belongs to the interior of this region (see Proposition 4-9 below). 

The set 3 ai in (4.20) contains the points x for which one pair of roots of P{z; a) coincides with the branch 
points ±io; these points are easily computed and are precisely the 2N -\- 1 points in (1.23). However we 
cannot positively conclude for general N that they are on the boundary of Pat, although this is quite evident 
from the numerics. Also, the detailed shape of dP^, beyond the easily established discrete Z 2 N +1 symmetry, 
is hard to describe in more detail; for example it is not obvious how to conclude that it consists of AN -\- 2 
smooth arcs for N G Z> 2 , as the Figure 1 clearly shows. We find it however already sufficiently interesting 
that we can narrow down the boundary of A as a subset of a simple set of equations (4.20), although we 
cannot completely describe it. 

4.4. At the center of Pn- In a small vicinity of x = 0, we have 

2N^l 

P4N+2iz)= Yl iz-al){z-af) 

k^l 

with for j ^ k. As in [5], the branch points a/" = —a)) are partially determined through (4.14),(4.15) 

and in addition through Boutroux type conditions 

g(z)dz]=0, 7 jePi(A,Z); y(z) = (P4 n-i-2 (z)) P 
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The latter are imposed on the hyper elliptic curve X = : nP' = P 4 Ar+ 2 (-z)} which is obtained by 

crosswise gluing together two copies of C\B with B = U^i ^ U^=i ® 2 /c+i]' 

Solvability of the resulting system for would now follow as in [5], but here we are only interested 

in the case x = 0. 


Proposition 4.9 (Compare Proposition 3.9 in [5]). For x = 0 the g~function is obtained from (4-13), (4-9) 
using 

= ^ 2 ]^ (^(^))' ’ = PiN+2iz) 

I _ 1 i-TTfc 

which is defined and analytic off z £ B with the branch points af = afc,o = 2 2 iv+i esfv+i , 2N + 1. 


= 2 


4N+2 


Proof. Local behavior. Near 2 = 0 we have 


viz) 


1 1 
2 02W+2 


(l + 0(z4^+2)) 


(4.21) 


and near infinity clearly y{z) = y + 0{z ^). Note that the determination of the root near z = 0 is the 
opposite. 

Boutroux condition. We have 



{R{z)) 


dz 

^2N-\-2 






-2V-1 


{R{ujz)) 


dz 

^ 2 V +2 


f'(^j+2,0 


{R{z)) 


dz 

^2N+2 ’ 


(4.22) 


and thus it is sufficient to verify the Boutroux condition 


for a specific j G 2A^}. This condition guarantees that all branch points lie in the zero level set of $. 

But for j = 1 it follows immediately that the integral is imaginary using the Schwartz symmetry. 

Connectedness of the level curves. First of all the set 4>(z) = 0 in C\{0} consists of one connected 
component alone; this is so because there are no saddle points and if there were two or more connected 
components, there would have to be a saddle point in the region bounded by them. We shall now verify 
that the level curves satisfy the necessary and sufficient conditions specified in Section 4.3. The critical 
trajectories must 

(1) connect all 47V + 2 branch points 

(2) obey the Z 4 Ar +2 symmetry because of obvious symmetry. 

A simple counting then shows that the only possibility is that exactly one trajectory from each branch point 
(in fact a straight segment) connects the branch points to 0 because the sign of $ changes 4A^ + 2 times 
around a small circle surrounding the origin. The other two trajectories must then connect the branch 
points. This is depicted in Figure 3. The discussion on the necessary and sufficient condition for the correct 
inequalities is now as explained in ([5], Section 3.1). □ 


Appendix A. Proof of the Miura relation (2.10) 

Remark A.l. We draw the reader’s attention to the various notations used in this section, 

t = (ti,12,13,^4, •■•) , to = (ti,0,t3,0,t5,...) , to = (ti,0, 2^73,0,2^^15,...), tjGC. (A.l) 

Let k,i G Z>o and introduce 

Bkit) = fz f z'^e’’^^) —, i7(z) = d(z;t) = ^n,e(t) = det [fie+j+k-2i't)]2k=v ^o.£(t) = 1 

771 Js z z 

where S C C denotes the unit circle traversed in counterclockwise orientation. Recalling (1.6) we see that 

= hk{t) 




4>1 
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Figure 3. The placement of the trajectories of the quadratic differential r] in the case of 
a: = 0; the red arcs are the branch cuts of the the g-function. The cases presented correspond 
to = 2, 3 (left/right). Also indicated is the foliation by the trajectories Cr of the region 
Koo (see Section 4.3). The shaded (cyan) regions indicate where $ > 0; all the trajectories 
that issue from the branch points constitute the set $ = 0. 


and thus with (1.5), 

A„,^(t) = S(£+„ + l)r.(t). 

In particular, by Lemma 3.1, we know that for the special value £ = 0 we have the identity 

A„,o(to) = S(„+i)>.(to) =2"”%5^(to). 

Next, let {pn,eiz)}n>o be the monic orthogonal polynomials associated with the measure 


1 Hr 

dMz) = —, zgS, (G Z>o; 

ZTTl Z “ 


Pn/(^)Pm,K^)dt'£(2) = KS„ 


, zGC\S; = -27ri 


A„ 7 (t) 


It is well known [10] that the matrix 

/ \ 771-1^ r ( \di'£,(w) 

satisfies a Riemann-Hilbert problem, i.e. r„ 7 (z) is analytic for z G C\S and we have the conditions 

Proposition A.2. The following identities hold for the Hankel determinants A„ 7 (t). 

A„ 7 +i(t —[z]) f-A A„ 7 _i(t + [z]) 

—-— = (-1) (r„7(z)j 


and 


A„ 7 (t) 

Ayi+i7(t) 

A„ 7 (t) 


A^2-1-17(1) 
le the 

z G S] r„ 7 (z) = + r„ 7 (oo)— + O (z 

aants A„ 7 (t). 

= (-i)"(r„7(z))2. 


(A.2) 


00. 


= -27ri (r/7(oo)) 


A„ 7 (t) 

A„i 7 (t) 


12 ’ 


A„7(t) 271 


= ^ (r/ 7 (oo)) 
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(A.3) 

(A.4) 


where [z] denotes the infinite vector of components (z, ’y, ^,...), i.e. 

f z^ z^ 

tT[z]= itiTz,t2 T yAsT y,--- 

Proof. The two identities in (A.4) follow simply by inspection of the expression (A.2). As for the identities 
(A.3), the proof follows from Heine’s formula for the orthogonal polynomials and the observation that 

exp [i9(w; t — [z])] = ii;^“^(w — z) exp [i?(w; t)] , exp [ 19 ( 11 ;; t + [z])] =-exp [i9(ii;; t)]. 
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Indeed, we have 


A„,£+i(t- [z]) = det [^ie+j+k-i{t - ^ ^ fl W)] ^ 

' S”’ • ^ j- ,•_1 


= ~ f = (- 1 )” det 

Kfc i=i 


i<fc 

Atf 


M^+n —1 

1 


i=i 

t^£+r, 


/^£+ 2 n-l 

Zn 


= {-l)"pn,e{z)An,l{t) 


where we used the well-known representation of orthogonal polynomials in terms of moment determinants 
(see, e.g. Proposition 3.8 in [10]). The second identity can be found in [4], but we can give here a direct 
derivation using Andreief’s identity [3]. Recall that Then 




2 dlZi{wj) 


j<k 


i=i 


= ® det 


..k-l 


= det 




,.fc-i 


J j,k=l 


det 




Wa 


n 

i=i 


dviiwj 


Wa 


J j,k-l 


det 


w 


k-l 


n n 


^Wj - 


Y[diye{wj 


(A.5) 


Multi-linearity allows us to replace the monic powers in the first determinant by the monic orthogonal 
polynomials pj/{w), so that we obtain 

..k — 1 -in n 


(A.5) = f det Pk-i^tiwi) 


det 






Wi — z 


-I 7=1 


diz^{wj). 


(A. 6 ) 


Now, in the second determinant we can subtract to the columns 2 < k < n the multiple z^ ^ /{wk — z) of 
the first column, thus obtaining 


(A. 6 ) = ® det Pj-i,i{wk) 


pk=l 


det 


Wi—Z Wi—Z 


Wn—Z W„—Z 


Using now Andreief’s identity we obtain 

/sPo,£(w)fEfdRKii;) ••• 
§sPdi{w)^,dvt{w) 

[A.l) = n \det 

_§sPn-lAw)^^i^ §sPn-lAA^A^M ■■■ 

but due to orthogonality the matrix above has the following structure 

§sPoAA^A^ % * ■■■ * 


ApoAA 

A Pi/A) 


Y\_dviAj) (A.7) 

3 = 1 


dvt{w) 

di3i{w) 


w—z 


§gPn-lAw)- 


-dvi{w)_ 


(A.7) = n! det 


LJS 


§sPiAA^A^ 0 /ii * -k 

: : 0 ■'. * 

: : : ■. /in -2 

p„_i.,(u;)^ 0 0 ••• 0 . 


= (-l)"+^n! d> pn-iAA 


di^iiw) 


w — z 


i =0 


However 


and therefore 


hj = ® pA{w)<^AA = 


Aj+i,^(t) 

Aj,«(t) 

di3i{w) 


An.^-2(t + [Z]) = (_l)»+l ^"-b^W I Pn-lAA^'A^ = (-l)"A„,,(t)(r„,,(^)) 

Ao,^(t) Js w-z 


22* 
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□ 


A.l. Dodgson-Hirota bilinear identity. Consider the following matrix-valued function 




gi5(z;t)-i?(z;s) q 
0 


r;;,+i(z;s), zec\(5u{0}). 


A direct inspection using the jumps of shows that this matrix has no jumps on the contour S', however 
an essential singularity at z = 0 due to the presence of the exponentials. We can thus compute the contour 
integral below in two ways. First by evaluation as a residue at infinity; 


27ri 


/ TT / \ dZ 

® HnAz-,t,s) - = 

'1 - (r7^-i(oo;t))i2(r7^+i('^;®))2i *’ 

/|z|=R ^ 

'k k 


(A.8) 


where the * indicates expressions which are not relevant to the steps below. Secondly we evaluate the left 
hand side in (A.8) as a residue at z = 0, but we are only interested in the (ll)-entry, 

i £l_^ t, ^ = res {ru,e-i{z; t)) (r„.f+i(z; s))^^. 

Hence with (A.8) and Proposition (A.2), 


res -e 

z=0 Z 


^ pi?(z;t)-i?(z;s) [^])An^^(s+ [z]) _ ^ ^ (t) (s) 


(l^)Ayj^£^_i(s) (t)A^^£^i(s) 

or equivalently 

res - [z])A„ ^(s + [z]) = A„,^_i(t)A„ £+i(s) - A„+i,f_i(t)A„_i ^+i(s). (A.9) 

z=0 Z 

Remark A.3. Identity (A.9) closely resembles a “Hirota” version of the classical Dodgson determinantal 
identity, for if we set t = s then (A.9) reduces to the Dodgson identity for Hankel determinants, 

^n,i = ^n,e-l ^n,l+l — ^n+l,e-l^n-l,e+l- (A.10) 

We now rewrite equation (A.9) with the substitution ti-At + h, s = t — h and define 

HDnj{t,h) = res A„_f(t + h - [z])A„_f(t - h + [z])^ 

— A„_£_i(t + h)A„_^+i(t — h) — A„+i^^_i(t + h)A„_i^^_|_i(t — h) 

so that (A.9) can be written in the compact form 

HDn^e{t,h) = 0 , Vt,h, (A.11) 

For the rest of this section we shall set all even times to zero, i.e. we choose t = to. Now use Corollary 3.2 
in conjunction with (A. 10), 

A^4(to) = A„_o(to)A„_2(to) — A„+i^o(to)A„_i_2(to) = 2(—l)"A„_o(to) A„+i_o(to), (A.12) 

and recall Lemma 3.1, 

A„,o(to) = S(«+i)"(to) = 2""%^^ (to). 

Hence with (2.9) for t 2 j+i = 0, j > N and ti = x, 

A„,o(to) = 2-"'e29A^;t)^ A„+i,o(to) = 2-^^+A" (A.13) 

Differentiating (A.11) with respect to hj we can derive a whole hierarchy of equations, however we are only 
interested in one particular identity: 


^HD„,£(to,h) 


h=0 


d^An+i^e-i . . An-i,E+i „ 5A„+i^^_i clA„_i_^+i 

— ^n+l,£-l - TTTTi -h Z 


dtl 


n,£-l . r, ^A„_f_i clA„_^+i 

■7^79 ^ P77 r 

utl uti 


dtl ' dtl 

d^An^£+i ^9^1nA„_^ 
dtj ^ dti^ 


dtl 

(A„,,)2 = 0 (A.14) 
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and the argument of all determinants in the right hand side equals t = tg. For £ = 1, with (A.14) and (3.4), 
(A.12) this leads to 

0 = (^ra+i,o^n,o + ^n+i,oA" g — 2A(j_,_j gg) + (In A„_|_i^g + In A„^o)" A„+i_oA„_g, (') = 


dti 


which can be rewritten as 

0 = 


A" A" A' A 

^ n +1,0 ^ n ,0 _ 2^11+1,0 ^ n ,0 

A„+i,o A„^o A„_|_i_g A 


— ^ + (In A„+i^g + In A„ g) 
n,0 / 


d , A 
A 


n,0 


n+1,0 


- 2^^2ln(A„,gA„+i,o)+ (^^^^In 

and after simplification with (A. 13), 

2dlgn{x]t) = -dlWnix;t) - (d^Wnix^t))'^ 
which completes the proof of (2.10). 


Appendix B. The outer corners of the regions Pn for N = 2,3 


In this section we offer a proof that the points (1.23) belong to the boundary of Pn- The proof is 
a verification that the inequalities for the effective potential are fulfilled at the particular values of a{x) 
determined in (4.8). These correspond in the a-plane to the points (1.23) in the cc-plane. The proof is a 
simple deformation argument starting from large |a| (and hence also large x). 

Observing various panes in Figure 1 and using the 'Z 2 N +1 symmetry of the region, it is sufficient to show 
that the point 


a 


[AT] 

0 


)(-!)» 



2iV + l 


AN] 


i-lf 


^( 2 A + 1 ) 


/2A + 1 

2N 




(B.l) 


(or rather its x-image) belongs to the boundary of Pn- This 
point is alternatively positive or negative, depending on the 
parity of N. Consider now in some detail the case N = 2; 
then ~ 0.944 (a;[f^ ~ 2.36021). In this case the polyno¬ 
mial P{z; a) (4.6) equals 


P = P2(z;a) =z'^+j^ 
4 


1 

2 a 


Let z^(a), j = 1,.. .2 denote the roots of P 2 - We know 
from the argument in Section 4.3 that for |a| large the in¬ 
equalities are fulfilled; as we deform a from larger absolute 
values to smaller ones, these inequalities can fail only if the 
sign of ^{z^{a);a) changes. 

We now simply have to verify that the sign of 4>(2;j(a); a) 
remains constant as a decreases from -too to the critical 

f2l 

value Oq (corresponding to x decreasing from -too to the 
f2l 

rightmost corner Xg ). Since the four roots admit an ex¬ 
plicit expression in terms of a, this verification is a simple 
exercise in calculus. To be more precise, one pair that we 

denote z^(a) is purely imaginary and lies on the zero level 

f2l 

set of $(z; a) identically for a € [og , 00 ); this is not a cause 
of concern because it belongs to the level curve (in fact a 



a 


Figure 4. The graph of the value of $ at the 
saddle point Zi- 


I fol 

straight line) joining 2 ; = ±ia to z = 0. The other pair zf^(a) is real for a G [flg , 00 ). Then one can easily 
verify that 


F(a) = 4>(zf(a);a) 
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is, depending on which of the two member of the pair, a monotone increasing/decreasing function of a € 

\ 2 ] 

[og , oo) and not changing sign. This verification uses Lemma 4.5 and the explicit expression for the roots, 
so that (for a real) 


d 

da 


$(zj(a);a) 


- 3 ? 




dx daP{z-, a) 
da P'{z\a) 



( \Jz^ ^ a? daP{z;a) 
y za P'{z-,a) 



(B.2) 


f2l 

In Figure 4 we display the graph of $(zi(a);a) in the range [cg ,oo); the monotonicity can be shown by 
inspecting the sign of (B.2); we leave the detail to the reader. The argument above can be repeated for 
N = 3, but for larger N we were not able to hnd a unifying argument. 
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